Low-Luminosity Extragalactic Water Masers toward M82, M51, 

and NGC4051 
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ABSTRACT 

Sub-arcsecond observations using the Very Large Array (VLA) are presented of low- luminosity 
H2O maser emission in M82, M51, and NGC4051. New maser features have been detected 
within the M82 starburst complex. They are largely associated with star-forming activity such 
as optically-identified starburst-driven winds, H 11 regions, or the early phase of star-formation in 
the galaxy. H 2 maser in M51 consists of blue- and red-shifted features relative to the systemic 
velocity of the galaxy. The red-shifted features are measured at the northwest of the nuclear 
radio source, while the location of the blue-shifted counterpart is displaced by ^2" from the 
radio source. A small velocity gradient closely aligned with the radio jet is detected from the 
red-shifted features. The red-shifted maser most likely amplifies the background radio continuum 
jet, while the blue-shifted counterpart marks off-nuclear star-formation in the galaxy. All of the 
detected maser features in the narrow-line Seyfert 1 galaxy NGC4051 remain unresolved by new 
VLA observations. Due to the low luminosity of the maser, the maser excitation is not directly 
related to the active galactic nucleus. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: individual: M82, M51, NGC4051 



1. INTRODUCTION 

Since the discovery of the first 22 GHz extra- 
galactic H2O maser towards M33 by Churchwcll 
et al. (1977), significant progress in studies of 
extra-galactic H2O masers has been made, with 
the first detection of a nuclear H 2 maser toward 
an active galactic nucleus (AGN) in NGC4945 
(|Dos Santos k Lepine Ifl979h . A number of single- 
dish surveys searching for new H2O masers in 
active galaxies have been conducted since, in 
which more than 1000 galaxies were observed. 
These surveys were stimulated by the VLBI imag- 
ing of H2O maser components in a sub-parsec- 
scale thin, warped, edge-on disk displaying Kep- 
lerian rotation around th e nucleus of the Seyfert 2 
galaxy NGC 4258 (e.g., iHerrnstein elaTl Il998l ). 
These surveys have increased the total number 
of extra-galactic EL O masers to ~80 at p resent 
(jHenkel et al. 11200a iKondratko et al.ll2006t ). 

It is generally recognized that extra-galactic 



H2O masers can be grouped acc ording; to their 
isotro pic luminosity (L iso ) (e.g., iGreenhill et al. 
19931 ). High-luminosity (L iso > 10 L Q ) H 2 masers 
are associated with active nuclei. Low-luminosity 
(L iso < 10 L Q ) H 2 masers are mostly associ- 
ated with star-forming activity however, some 
of them may also contain or be in low radio 
luminosity AGNs. The luminosities of Galactic 
H2O masers that are commonly observed in the 
envelopes of evolved stars and star-forming re- 
gions range, typically, from 0.001 L to 1 L Q , 
and no H2O maser with a luminosity well above 
1 L Q has ever been found our galaxy. Ac- 
cordingly, H 2 masers observed in nearby star- 
forming or starburst galaxies such as M33, M82, 
NGC253, and NGC6946 with Upvalues of 0.01- 
1 Lq ( e.R., Churchwell et al. 1977t Claussen et al 



1983 iHo et all Il987l : iBaudrv fc Brouilletl 11996) 



are most likely to originate in star-forming ac- 
tivity in the host galaxy. On the other hand, 
there has been debate over the nature of a sub- 
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population of H2O masers with luminosities of 
1 L0 < Lj so < 10 L©, most of which are consid- 
ered to arise in prominent sites of star-formation 
or starburst-activity in galaxies. However, some 
of their host galaxies show AGN activity; low- 
luminosity water masers have been detected to- 
wards LINER or type 1 Sey fert nuclei in recent 
sensitive single-dish s urveys ( Henkel et al.l 120021 : 
Hagiwara etaD 120031 : iBraatz et al.l l2004h . They 
could simply be low-luminosity analogues of the 
high-luminosity masers in narrow-line AGNs such 
as LINERs and type 2 Seyferts, and lower maser 
luminosities might be explained by a close to face- 
on view of obscuring structure around the nucleus, 
such as a disk or torus in the line of sight, or a 
misalignment between the disk and a nuclear con - 
tinuum in the line of sight (jHagiwara et al. I l2003h . 
It is important to explore the sub-population of 
extragalactic masers with 1 L Q < Lj SO <10 L Q at 
high angular resolution, as they could advance the 
study of extra-galactic star- formation or probe the 
nuclei in low radio luminosity AGNs. 



is 



M82 (D = 3.63Mpc: lFreedman et al 
a well-known nearby starburst galaxy with a long 
observing history at various wavelengths. The 
galaxy hosts a number of exotic radio sources 
within its starburst complex, including main- and 
satellite-line OH masers and absorption, which 
trace cold molecular material in the galaxy (e.g., 
Seaquist et "aTlll997t [McDonald et alJl2002n . Sev- 
eral low-luminosity H2O maser features were de- 
tected at 1.4" resolution in earlier Very Large 
Array (VLA) observations (jBaudrv fc Brouilled 
1996j) however, the resolution was not sufficient to 
pin down or resolve the maser to investigate its 
association with radio continuum sources such as 
the nucleus, jet, compact Hn regions, or super- 
nova remnants (SNRs). 

The Whirlpool gal a xy, M 51 (D = 9.6 Mpc; 
Sandage fc Tammann ( 19751 )). has provided an 
ideal laboratory for studying molecular materi- 
als in a galactic star-forming environment. These 
materials are more abundant in the spiral arms 
than in the central region. The nuclear re- 
gion is very complex, showing molecular emission 
that is asymmetric in both position and velocity 
with respect to the nucleus: the dominant red- 
shifted emission peak is 1 " to the west of the ra- 
dio/optical nucleus (e.g., Scoville fc Youngj 19831: 
Scoville et aHll998l ). Ho et al. (1987) first reported 



the discovery of low-luminosity H2O maser emis- 
sion towards the galaxy. A snapshot observation 
pin-pointed the location of the known red-shifted 
maser features (jHagiwara et al.l l200ll) northwest 
of the continuum peak of the galaxy, while the 
blue-shifted feature(s) have not previously been 
measured. Due to the insufficient velocity cov- 
erage of the earlier observation, the overall dis- 
tribution of the maser emission has not yet been 
explored in detail. 

[Adams I 



NGC4051(D = 9.7 Mpc; lAdams I (Il977l) ) 
is known to be a narrow-line Seyfert 1 galaxy 
(NLS1) showing 'broad' emission line in H/3 
(FWHM < 2000 km s _1 ), narrower than that typ- 
ically observed (FWHM = 2000-10000 km s" 1 ) to- 



wards broad-line Seyfert 1 galaxies (jOsterbrock fc Pogge 
Il985h . NLSls are characterized by their housing of 
low-mass black holes radiating near t heir Edding- 
ton limits (e.g.. IWilliams et al.ll200l . NGC4051 
exhibits very strong X-ray intensi ty variability on 



many different timescales (e.g., lLawrence et al 



19871 ). which is interpreted as a result of its nu- 



clear activity and extended starbur st component s 
dominated by s o ft X-r ay emission (jSingh lll999h . 
Hagiwara et al. (|2003l) detected an H2O maser 



in NGC4051 and identified it towards the cen- 
ter of the gala xy with a complex nuclea r radio 
structure (e.g., lUlvestad fc Wilson! 1 19841 ). The 
maser luminosity is estimated to be ~ 2 L Q . The 
maser features span ^300 km s _1 and straddle 
the systemic velocity of Vlsr= 730 ± 3 km 
s (jHagiwara et al. 2003) nearly symmetrically. 
These Doppler-shifted maser features appear to 
indicate the presence of a rotating disk however, 
there has been no compelling ev idence found to 



supp ort the presence of a disk (jHagiwara et al 



2003) . The origin of low- luminosity H2O masers in 



NLSls, including NGC4051, remains unexplored. 

In this paper, I report sub-arcsecond imaging 
of H 2 masers towards these three nearby galax- 
ies, in order to update the earlier observations 
with higher resolution and sensitivity and to fur- 
ther explore the nature of low-luminosity masers. 
All of them are so weak (flux density < 100 mJy) 
that they cannot be studied in VLBI observations. 
Thus, sub-arcsecond imaging is currently the best 
method to resolve the masers. In section 5 the 
statistical infra-red and radio properties of the 
H 2 masers are discussed. Throughout this pa- 
per H = 75 km s _1 Mpc _1 is adopted. 
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2. OBSERVATIONS 

Spectral line observations at 22 GHz using the 
NRAC0 VLA were carried out for measurement of 
H2O maser emission (6i6-523 transition) towards 
three nearby galaxies M82, M 51, and NGC4051 
from 2002—2003. All the observations were per- 
formed on these galaxies when the array was in 
the A Configuration. All the observations were 
made by employing two intermediate frequency 
(IF) bands of width 12.5 or 6.25 MHz with a single 
polarization, divided into 32 or 64 spectral chan- 
nels, yielding velocity resolutions of 1.32 or 2.63 
km s -1 per channel. Each IF band was centered 
on a velocity selected to cover all the maser emis- 
sion, with respect to the local standard of rest 
(LSR). 

Because low-luminosity extragalactic maser 
emission is difficult to detect within an at- 
mospheric coherence time at 22 GHz, phase- 
referencing observations were employed, using a 
nearby calibrator source (Table 1) several de- 
grees away from each galaxy. The observations 
were thus performed by switching to the phase- 
referencing source, typically every 2 minutes for 
30 s. Amplitude and bandpass calibration were 
performed using observations of 3C286, and the 
flux density scale was estimated to be accurate 
to within 10%. The data were calibrated and 
mapped in the standard way using the NRAO As- 
tronomical Image Processing Software (AIPS). A 
summary of the observations, including IF veloc- 
ity ranges, synthesized beam dimensions, and rms 
noise for each observing run, is provided in Ta- 
ble 1. After the phase and amplitude calibrations, 
the image cube was continuum subtracted using 
channels with no significant line emission in order 
to obtain only the line-emission. Thus, the maser 
emission in M 82 was separated out from the con- 
tinuum emission. The continuum image of M 82 
was produced, collecting only line-free channels 
within a single 6.25 MHz IF band. In the (u, v) 
data of M51 and NGC4051, no continuum emis- 
sion was detected to a 3 a rms noise level of 0.6 
mJy beam -1 for M51, and 0.48 mJy beam -1 for 
NGC4051. 

In this paper, the B1950.0 coordinate system is 



adopted for the observations and figures of M82, 
since this system was used in many earlier works 
on M82 in the literature. 

3. RESULTS 
3.1. M82 

Figure [T] shows the locations of the detected 
H2O masers, and other known sources in M82. 
A single-dish spectrum and individual VLA spec- 
tra of the H2O masers are displayed in Figures 2 
and 3. All five VLA spectra had detections above 
the ^5 er level. However, one might suggest that 
the detected features at N2 and N3 need to be 
confirmed. All of these maser features remained 
unresolved at the angular resolution of ~93mas, 
corresponding to 1.7 pc at the adopted distance, 
in the robust weighted maps. 

The entire galaxy was searched, and no other 
water maser emission was found in the velocity 
range of Vlsr= 50—150 km s -1 . The known red- 
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shifted emission (e.g. JClaussen et al.lll984h in the 
galaxy was not seen in this observation. In addi- 
tion to the four masers S1-S4 (us ing the designa- 
tion by IBaudrv" fc Brouilletl (|l996h . detected with 
the VLA in the C-configur ation at 1.4" resolution 
(|Baudrv fe Brouilletl Il996h . mree maser sources 
are newly detected in our observing run (Fig. [3]). 
They are hereafter labeled Nl, N2, and N3, respec- 
tively. None of the positions of the H2O masers 
coincide with those of the OH masers, luminous 
X-ray sources, and infrared peaks that have been 
measured to date. 

Table [2] shows the positions, peak flux in- 
tensities, velocities, and velocity widths of the 
detected masers at the highest angular resolu- 
tion of 0.93". These new maser positions are 
constrained better than those at 1.4" resolution 
me asured with the VLA in t he C configuration 
by iBaudrv fc Brouilletl (|l996h . The shapes of 
th e detected maser spectra that were reported 
in IBaudrv fc Brouilletl (|l996l ) were relatively un- 
changed since their observation in 1993. The 
positions of SI and S2 are consistent to within 
0.08" between the two observations, and the dif- 
ference is most likely due to different interfer- 
ometer beam sizes. No maser emission at or 
near the location of S3 and S4 was detected to 
a 3a level of 3.3 mJy beam -1 in this observa- 
tion, probably due to intensity variation of the 
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maser. Table [2] lists the positions of the maser 
features at S2, in which the emission profile is 
peaked at several velocities, spanning ss 30 km 
s -1 , but all lying in an ~ 0.01" region cen- 
tered on R.A. (B1950)=09 h 51 m 42. 1812 s , decl. 
(B1950)=69°54'59.260". 

Table [3] lists the nearest discrete continuum 
source and radio recombination line (RRL) for 
each H2O maser. In order to compare the posi- 
tions between the masers and discrete cont i nuum 
sources in the literature (McDonal d et al. 2002, 



0.01", centered on R.A. (J2000) = 13 h 29 m 52 s .709 
decl. (J2000)=+47°ll'42".79, wh ich is consistent 



and references therein), the published positions 
of two known sources at 15 GHz were adjusted to 
those on our 22 GHz continuum map, which led 
to a sensible shift of < 0.05" for 41.95+57.5 and 
43.30+59.2. Therefore, one can estimate that po- 
sitions between the two observations can be com- 
pared at an accuracy of < 0.05" . The maser posi- 
tions in Table [3] show an offset of 0.2-0.3" (corre- 
sponding to 0.36-0.54 pc) from the nearest contin- 
uum peaks. Such a large separation makes it diffi- 
cult to conclude that the detected maser is associ- 
ated with an Hll region or SNR. Thus, the masers 
are not essentially associated with any continuum 
source at this resolution. Alternatively, (compact) 
H II regions were not detected due to the sensitiv- 
ity limit of these observations. 

3.2. M51 

Figure @] displays VLA spectra of H2O masers 
detected towards the central region of M51. The 
red-shifted features, spanning Vlsr= 538 - 586 
km s , were prev iously detected i n VL A snap- 



shot observations (jHagiwara et al.l 120011 ). One 
of the known blue-shifted features, peaking at 
Vlsr = 445 km s _1 , is detected by interferom- 
etry for the first time. These Doppler-shifted 
velocity features straddle the systemic veloc 



Scoville fc Young 



i ty (V sys ,LSR = 472+3 km s" 
(|l983l) 1 of the galaxy. 

Figure [5] shows the approximate locations of 
the blue- and re d-shifted features on the 8.4 GH z 
VLA-A image of lBradlev. Kaiser, fc Baan (2004). 
None of these features is resolved by the uni- 
formly weighted beam of 82 x 77 mas 2 . Assum- 
ing the masers have Gaussian distributions, the 
beam-deconvolved size of the red-shifted maser 
is ^25 mas, corresponding to 1.2 pc at the ac- 
cepted distance to the galaxy. The positions of 
all the detected red-shifted features are within 



with the earlier observation by lHagiwara et al 
(120011) . The blue-shifted feature is located at R.A 



(J2000.0)=13 h 29 m 52 s .549, decl. (J2000.0)=+47°ll'43".34. 
The parameters of the maser features are listed in 
Table H 

The projected distance between the blue- and 
red-shifted features corresponds to 82 pc (or 1.7"); 
the blue-shifted feature is well separated from the 
other known features and so they are not kine- 
matically related. On the other hand, in both this 
and the earlier VLA observation, the red-shifted 
features are observed ~180mas (or 8.6 pc) north- 
west of the radio con tinuum peak of the galaxy 
(jHagiwara et al.ll200ll ). suggesting that the maser 
does arise from the radio jet rather than in a low- 
luminosity nucleus of the galaxy. 

Figure [5] also displays the first moment map 
produced from the spectral-line data cube con- 
taining only 10 velocity channels of red-shifted 
features in Fig |4j A weak velocity gradient, 
roughly from south to north is seen along the axis 
PA. = 155°. The value of the gradient is, approx- 
imately, 10km s _1 per 0.1", or 2km s _1 pc -1 . 
The gradient is very weak but real, at least from 
554 to 563 km s _1 , although the two extrema at 
the edges Vlsr = 564 and 553 km s _1 are less cer- 
tain. However, the zeroth-moment map is well 
sampled with a good signal-to-noise ratio toward 
the center portions of the source. As the VLA 
beam is 82 x 77 mas at P. A. = 6°, there is no pos- 
sibility that the gradient along P.A. = 155° is aris- 
ing from effects of the beam shape or orientation. 
Accordingly, I conclude that the observed velocity 
gradient and its axis are real. 

3.3. NGC4051 

VLA spectra of the H2O maser towards the 
center of NCC4051 (Fig. ^ cover most of the 
known features in the range Vlsr = 640-780 km 
s _1 . The location of the systemic feature peak- 
ing at Vlsr = 717.6 km s _1 is measured at R.A. 
(J2000) = 12 h 03 m 09 s .610, decl. (J2000)=+44°31'52".68, 
and all other features are confined to within about 
0.01" (or 0.5 pc) from the systemic feature, result- 
ing in the masers being unresolved at this reso- 
lution. The beam-deconvolved size of each maser 
feature is smaller than 40 mas, that is, about 2pc. 
The result is consistent with the earlier VLA-A 
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snapshot by IHagiwara et al. ( 2003f h Although 
this new observation has improved (u,v) coverage 
over the previous single VLA snapshot, no new 
useful information is obtained. 

4. Discussion 
4.1. M82 

4.1.1. H% Masers and Luminous X-Ray Sources 

The positions of the luminous hard X-ray 
sources in M 82 were measured within an ac- 



curacy of 0.7" by Ch andra ( Kaaret et al.l 12001 



Matsumoto et al.l 12001 ): however, none coincide 
with those of the H2O masers detected in this ob- 
servation. In addition, there was no high bright- 
ness radio source detected at 22 GHz towards 
any of these luminous X-ray sources, suggesting 
that these X-ray sources are not likely to be low- 
luminosity AGNs. Since the masers do not overlap 
any high-brightness radio sources like an AGN or 
jet, they are off-nuclear masers that are not di- 
rectly amplifying the background continuum, as 
in the case of nuclear H2O masers. This is con- 
sistent with the fact that the maser luminosity in 
the galaxy is 1-5 orders of magnitude lower than 
those of the nuclear H2O masers. 

4-1.2. H2O Masers and Other Sources 

None of the H2O masers detected in this 
observation are associated with infra-red peaks 
at 2.2 /im, which were considered to be low- 
luminosity AGN ca ndidates or the dy namical cen- 
ter of the galaxy ( Dietz et al. 19861 ). suggesting 
that the masers in M82 are not nuclear masers. 
The separations between each H2O maser peak 
and the center of the nearest continuum source, 
given in Table 3, indicate that there are no 
H2O masers at the locations of the Hil regions 
or radio SNRs. According to the 5 GHz MERLIN 
and 15 GHz VLA-PT ob servations of continuum 
sources in the galaxy by iMcDonald et al I (|2002l ). 
the angular size of each Hn region appearing in 
Table 3 is - 0.1"-0.2"(2-4pc) in radius. It is un- 
derstood that the sizes of Galactic compact Hn 
regions excited by central stars are typically 0.05- 
0.5 pc and Galac tic H II regions are less th an ~ 5 
pc in radius (e.g. JGaray and Lizano lll999l ). Given 
the observed separations between the masers and 
the centers of compact Hn regions that range 



from 0.17" to 0.51", corresponding to 3-10 pc, it 
is less plausible that the masers in the galaxy are 
directly associated with Hn regions or ionizing 
central stars in the galaxy. In contrast, VLBI 
observations revealed that the H2O maser in the 
nearby star-forming g alaxy M33 is ass ociate d with 
the HII region IC133 (jGreenhill et al.lll993l ). It is 
thus plausible that some other H2O features in 
M82 might be found in Hn regions. The de- 
convolved angular sizes of RRL emission (H92a) 
in the galaxy displayed in Table 3 range from 
1.4"-1.7"(25-30pc), just greater than the VLA 
synthesized beam of 0.9" ((Rodriguez et al.l l2004h 
that includes the masers at SI, S2, and N2. The 
H2O maser would not have originated from the 
RRL -emitting ionized gas, but the RRL marks 
the prominent star-forming regions in the galaxy, 
where maser excitation actively occurs. 

The velocity deviation of the ionized gas (100 
- 130 km s _1 ) and molecular gas (110 - 145 km 
s _1 ) from the galaxy's rigid rotation velocity 
has been observed, and accounts for the expan- 
sion of starburst-driven winds in the galaxy (e.g., 
iMatsushita et al1l2005l ). but the kinematics con- 
necting the masering clouds with the ionizing 
medium is not obviously seen from this obser- 
vation. 

4.I.3. What Are the Masers in the Galaxy? 

The interpretation of the nature of the H2O masers 
in M82 is not straightforward, since none have dis- 
tinct co ntinuum counterpa rts as effective pumping 



agents. Wei fi et al ] (Il999l) imaged an expanding 



shell-like structure (superbubble) traced by CO 
and CO isotopes in the galaxy, with a diame- 
ter of !=al30pc and expansion velocity of ~ 45 
km s _1 . According to their analysis, the veloc- 
ity of the approaching side of the expanding shell 
is Vlsr= 95 km s _1 and the receding velocity 
is Vlsr= 190 km s _1 . Thus, the former has an 
H 2 maser counterpart at Vlsr= 86 and 94 km 
s _1 , but the latter is not covered in this obser- 
vation. If the maser at S2 is associated with the 
shell, it could be excited at the shock front of the 
starburst-driven winds. It is interesting to note 
that thermal 100 GHz continuum emission has 
been imaged inside the CO superbubble, and the 
distribution of the H2O masers roughly underlines 
the conti nuum emission originati ng from starburst 
activity (jMatsushita et al.ll2005h . About 10 OH 



5 



main-line masers are observed that are associated 
with b lue-shifted parts o f expanding wind-driven 
shells dWeifi et al.lllQQflh . and the velocities and 
locations o f these OH and H2O masers are largel y 
consistent (JPedlar Il200l iMatsushita et al.ll2005l ). 
These OH masers trace the 1667 MHz OH ab- 
sorption along the major axis of the galaxy, and 
the OH dynamics understood as soli d-body rota- 
tion a grees well with the CO emission (|Shen k. Lo 



disks around extra-galacti c YSOs that have n ot 
been studied to date (e.g., Torrelles et al. 1998f ). 



19951) and Hi absorption (jWills et al.ll2000f) . How 



ever, there are disti nct deviations from this at 350 
and - 70 km s" 1 (|Pedlar et al. 1 12004 ) . although 
neither is observed in our observations. Therefore, 
the detected blue-shifted H 2 masers are concen- 
trated on the same molecular disk that is traced 
by OH absorption and CO emission; however, the 
association of the H2O masers with the expanding 
shells or molecular outflows is not obvious from 
our data. 

Similarly, in the starburst galaxy NGC 253, the 
locations of low-luminosity H2O maser (~ 1 Lq) 
are measured towards th e nuclear region wi th the 
VLA at ~ 1" resolution (jHenkel et al.ll2004) . The 
association of the strongest blue-shifted maser fea- 
tures with the RRL, an expanding supernova, or 
starburst-driven winds is proposed , like the case 
in M82; however, there is no compelling evidence 
for that from their observations. 

The extent of the CO(2-l) molecular clouds 
and outflows is not inconsisten t with that of the 
H 2 maser (jWeifi et al.ll200ll) . Molecular out- 
flows in our Galaxy are observed in H 2 emis- 
sion, typically on scales of 1-10 AU, which is at 
least 100 times smaller than the masers measured 
on the parsec-scale with the VLA. Accordingly, 
even if the masers probed part of the outflow, they 
would not be resolved on these scales by analogy 
with the cases of young stellar objects (YSOs) in 
our Galaxy. It is believed that H 2 masers indi- 
cate sites of star-formation and appear at a cer- 
tain stage of the evolutionary history of proto- 
stars. It is also understood that CH3OH, OH, and 
H 2 masers are observed at different evolutionary 
stages of star formation, and H 2 emerges at the 
earliest stage of massi ve star formation d uring the 
rapid accretion phase ( Churchwell 2002). Masers 
in a galaxy that do not accompany compact Hll 
regions and OH masers may be the signposts of the 
early stage of star formation and are most likely to 
be associated with molecular outflows or accretion 



4.2. M51 

4-2.1. Blue-shifted Features 

The precise location of the known weak and 
variable blue-shifted features lyi ng from Vlsr 



435-4 45 km s in this galaxy (jHagiwara et al 



20011 ) has been an open question. One of them has 
been pinpointed in this observation for the first 
time. Taking the fact that the maser luminos- 
ity of the blue-shifted emission is approximately 
0.17 Lq from the VLA spectrum and is not associ- 
ated with any known radio continuum sources in 
the galaxy, one can infer that AGN activity is not 
the major source giving rise to the maser. The 
blue-shifted emission is displaced 1.7", or 82 pc, 
from the red-shifted counterpart near the nucleus, 
thus, no hint of physical connection exists between 
these velocity clusters. From what do the features 
arise? The peak velocity of Vlsr= 445km s _1 is 
i n the velocity rang e of the central HCN emission 
( Kohno et al.lll996l) . The features lie in the cen- 



tral dense region of the HCN emission, where the 
mean velocity field of the HCN emission is Vlsr~ 
480 km s" 1 , that is , 35 km s _1 shifted from the 
maser ( Kohno et al. 1996). Thus, it is impossi- 



ble to find common mechanics connecting the two 
different velocity gradients traced by the maser 
and the HCN. All that one can speculate is that 
the maser traces the most dense part (> 10 7 cm~ 3 ) 
of the circumnuclear molecular gas within 10 pc 
of the nucleus. As a result of this observation, 



the bi polar-jet model proposed in lHagiwara et al 



(|200ll ), in which the blue-shifted features are as- 
sociated with the approaching side of the jet, is no 
longer eligible. 

4-2.2. Central Kinematics 

One of the most exciting results from this obser- 
vation is that part of the velocity structure of the 
dense molecular gas within a few parsecs of the nu- 
cleus is resolved at sub-arcsecond resolution. The 
detected velocity gradient is about 2 km s" 1 pc^ 1 , 
and the direction of the gradient (PA ~ 155°) is 
similar to that of the radio jet but not to the ro- 
tational axis of the inner torus with a radius of 
70 pc (PA - 160°-165°) postulated from the dis- 
tribution of the HCN (J = 1—0) emission observed 
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at 4" resolution (jKohno et al.lll996l ). Most of the 
red-shifted maser features in Fig. |4] are neither in 
the velocity range of CO (J = 1-0) (Vlsr= 387- 
547 km s- 1 ) nor in that of HCN (J = 1-0) (V LS r= 
377-546 km s~ 1 ) . The velocity ranges of the maser 
do not agree with those of the HCN emission and 
the H CN torus model proposed in iKohno et al.1 
(jl996h . however, the axis of the HCN torus is not 
inconsistent with the direction of the maser veloc- 
ity field. 

Given the fact that the blue-shifted features are 
separated by ~80pc from the rest of the features, 
and do not coincide with the nucleus, the whole 
maser system is unlikely to trace a subparsec- 
scale Keplerian disk observed with the nuclear 
H2O maser, if it exists. However, it is possible that 
only the red-shifted features trace a part of a rotat- 
ing disk, and the blue-shifted features on the disk 
are invisible. Assuming that the observed maser is 
on a disk at a radius r from the nucleus and with a 
rotating velocity of V ro t , the observed velocity gra- 
dient (dV/dl) along P. A. = 155° can be expressed 
as dV/dl ~ d(V ro t 8 sim)/d(r6>) = V rot /r, where 
I is the projected distance, i is the disk inclina- 
tion, and the approximation of sin6> w 9 is made. 
Adopting the disk inclination of 80° from the edge- 
on cir cumnuclear- nuclear disk with a radius of 
70 pc (jKohno et al.lll996h . together with the ob- 
served parameters of dV/dl = 2 km s _1 pc -1 and 
V re d — V sys — V rot — 100 km s _1 , r is estimated 
to be 50 pc. The result suggests that the observed 
maser probes a disk on a scale of 10 pc but not 
a thin disk wi th r <0.1pc, such as the one in 
NGC4258 (e.g. lHerrnstein et al.lfl998h . 

Does the maser disk lie in an inner part of the 
circumnuclear HCN disk? Taking this value for 
the radius of a molecular disk, the mass confined 
within the disk is 1.2 x 10 8 M®, that is, 1 or 2 mag 
larger than the molecular hydrogen and dynami- 
cal mass within r < 70 pc, based on the estim ation 
from the HCN intensities (jKohno et al.lll996f) . Ac- 
cordingly, we cannot account for the detected ve- 
locity field with a subparsec-scale masering disk or 
a larger scale molecular disk. 

It is interesting to note that the velocity gra- 
dient of 0.65 km s _1 pc -1 along the radio jet 
axis was detected in CO(J = 3-2) at 4"resolution 
160 pc) using the S ubmillimeter Array by 



than that in this VLA observation. The CO ve- 
locity gradient in their position-velocity map cov- 
ers that of the maser at Vlsr= 555-565 km s _1 . 
However the observed CO gradient is more domi- 
nant in the direction perpendicular to the jet, and 
such a small velocity gradient in the thermally 
excited molecular gas is likely d ue to the internal 
turbu lence, as also mentioned in lMatsushita et all 
((2004J). Thus, there is no compelling evidence 
that the maser and CO(J = 3-2) trace the same 
kinematics, but both of them could be distributed 
in the same molecular cloud. In the VLBI ob- 
servation of the 1667 MHz OH maser towards the 
type 1 Seyfert nucleus in Mrk231, a weak veloc- 
ity gradient of 1.4km s _1 pc -1 across the torus 
structure with a radius of 65 pc is reported in OH 



emission (jKlockner et al.l j2003) , which does not 



agree with any thermal molecular gas structure. 
In this sense, the OH maser in Mrk231 is similar 
to the H2O maser in M51. The axis of the OH 
torus in Mrk231 is not aligned with the major axis 
of the nuclear continuum source; rather, it is per- 
pendicular to the nuclear continuum axis. This is 
different from the case of M51, in which the veloc- 
ity field of the maser is almost aligned with the jet 
axis, which might rule out the presence of a maser- 
ing disk in the galaxy. Note that there are still 
some other blue-shifted features whose positions 
have not yet been pinned down, some of which 
could be on a masering torus within 1 pc from 
the center, along with red-shifted counterparts. If 
the maser in M51 is resolved down to scales of 
0.1 pc using VLBI, the velocity structures of the 
maser will be interpreted differently. The most 
plausible explanation for the maser and its veloc- 
ity gradient based on these observations is that 
the maser arises in a dense molecular environment 
in the foreground of a radio jet, amplifying the 



backg round radio continuum jet (jClaussen et al 



1998), or, alternatively is from a s hocked dense re 



gion near the boundary of the jet (jHagiwara et al 
200 lj ). The weak radio intensity of the jet would 



Matsushita et al. ( 2004f ). although the sampled 



area in their image is more than 100 times larger 



account for the low luminosity of the maser. 

4.3. NGC4051 

The luminosity of the H 2 maser in NGC4051 
is ~1Lq, very low for nuclear masers. It 
might appear that such a low-luminosity maser 
is no t related to dir e ct A GN activity. How- 
ever, iHagiwara et alJ (|2003l ) proposed that the 
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low-luminosity of the maser and narrower to- 
tal velocity span of the Doppler-shifted fea- 
tures (~100 km s _1 ) might be due to low-gain 
maser amplification, resulting from a small in- 
clinati on of a disk-like struct ure surrounding the 
AGN. Hagiwara et al.l ( 20031 ) hypothesized that 
the maser originated from a less edge-on disk-like 
configuration surrounding a nucleus, which would 
cause relatively short gain paths in the line of 
sight. This cou ld account for t h e low -luminosity 
of the maser. iMadeiski et al. (2006), however, 



argue that the maser in NGC4051 is unlikely to 
be connected with direct AGN activity, in that 
the X-ray absorbing medium is significantly ion- 
ized in NGC4051, and such an ionized absorb- 
ing medium is unlikely to be physically close to 
the masering medium. They suggest that the 
origin of the maser in NGC4051 is in nuclear 
winds, because the broadly spread (~ 280 km s _1 ) 
narrow -line spectrum is similar to that of the 
Circinus galaxy. Circinus hosts highly Doppler- 
shifted maser features with a total velocity range 
of -160 km s" 1 that are introduced as nuclear 
wind components, together with features tracing 
an edge-on disk wi th a sub-Keplerian rotation 
(jGreenill et al.l I2003T ) ■ However, a physical mech- 
anism driving the nuclear winds in the galaxy has 
not yet been proposed. 

One can speculate that the medium giving rise 
to the maser is geometrically separated from the 
X-ray absorber. The inner edge of the absorber 
could be illuminated by direct X-ray radiation and 
hence ionized, which would also cause a high time 
variability of the X-ray-measured column density 
Nh, while the masering medium could exist in 
the outer edge of a thin disk that is much less 
ionized because of a large radial distance from the 
X-ray heating source. The known periodic inten- 
sity variability at hard X-ray ban ds might occur in 



the i nner edge of the disk (e.g.. iHerrnstein et al 



1998]), while the time-interval of the intensity vari- 
ability of the maser is very different. In the case of 
NGC 4388, the radius of the absorbing medium is 
estimated from the variability of Njj to be a few 
hund red Schwarzschild radii (R sc h) (|Elvis et al.1 
2004), that is 100 times smaller than the (sub- 
)parsec-scale obscuring medium observed in the 
H2O maser. H2O masers in NGC 4258 are dis- 
tributed in a region between 40,000 and 80,000 



R sr h from the nucleus (e.g., IHerrnstein et al 
ll998l ). Therefore, it is less plausible that the 
masering medium and X-ray-ionizing structure 
occur in a single physical structure. Thus, we still 
cannot rule out the presence of a maser disk in 
N GC 4051. In any case, the hypothesis proposed 
bv lHagiwara et al. ( 20031 ) should be tested at mil- 
liarcsecond angular resolutions by resolving the 
maser distribution on scales of 0.1 pc or less. 

5. FAR-INFRARED VS RADIO INTEN- 
SITY CORRELATION DIAGRAM 

There is a known strong linear correlation be- 
tween the far-infrared (FIR) flux density and the 
radio flux densi ty at 1.4 GHz from nuclear star- 
burst galaxies ( Helou et al.1 Il985 ; Condon et al 



19911 ). which implies that FIR radiation and non- 
thermal radio emission at 1.4 GHz are closely 
connected with star-forming activity. In Fig [7] 
there are three different comparisons: the ra- 
tios of the 1.4 GHz radio fluxes to FIR flux from 
high-luminosity H2O masers, as well as from 
low-luminosity H2O masers and prototypical OH 
megamasers (OHMMs) from ultra-luminous FIR 
galaxies. The low-luminosity masers and OHMM 
samples follow well the known FIR-radio corre- 
lation. The mean constant ratios of these three 
samples are also displayed in Fig [7j from which 
one can see that the ratio for the low-luminosity 
masers is steeper than those of the other two. 
The reason for this deviation is that the FIR 
fluxes are relatively more dominant for the low- 
luminosity masers, which supports the idea that 
the low-luminosity masers are powered by star- 
forming activity rather than AGN activity. We 
cannot discriminate the synchrotron emission at 
1.4 GHz from cither SNR or AGN ejecta, but the 
FIR emission is characterized as re-emission from 
the dust heated by starburst /star- forming activ- 
ity (e.g-. lHagiwara et al.ll2002l ). It is interesting to 
note that there is a population of high-luminosity 
masers that shows a higher ratio of the radio to 
FIR fluxes in the bottom phase-space in the di- 
agram, which may suggest a close connection of 
maser excitation to strong radio flux from AGN 
activity. 
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6. SUMMARY 

Low-luminosity extragalactic H 2 masers 
have been identified in different physical envi- 
ronments related to star- formation, such as Hll 
regions, YSOs, or possibly SNRs. Given the 
fact that the luminosity of the strongest Galac- 
tic H2O maser, W49N, is ~1Lq, the origin of 
these low-luminosity masers can be explained by 
star-forming activity in their host galaxies. Many 
of the host galaxies of these masers contain star- 
forming regions or exhibit starburst activity but 
do not contain an AGN. These H 2 masers may 
reveal new stellar phenomena relevant to extra- 
galactic star-forming activity. 

The kinematics of H2O masers in M82 is 
broadly consistent with OH or CO solid body 
rotation along the galaxy's major axis, but lo- 
cally the masers are not associated with any other 
molecular emission. Including the three new de- 
tections, the masers are not directly associated 
with any continuum sources , such as compact H 11 
regions or radio SNRs, suggesting that the masers 
arise from the earliest stage of star-forming ac- 
tivity in the galaxy. For the masers in M51, the 
detected blue-shifted features are significantly off- 
set from the nuclear radio continuum, while the 
red-shifted counterparts are in the nuclear radio 
source. Seemingly, the maser in the galaxy am- 
plifies the radio jet continuum, which needs to be 
confirmed at higher resolution. The most remark- 
able result is that a small velocity gradient roughly 
along the jet has been detected from red-shifted 
emission at 0.1" resolution. However, the inter- 
pretation of the velocity gradient is not straight- 
forward. The nature of the maser in NGC4051 
is controversial and also needs to be studied at 
higher resolutions. A study of the distribution 
of the rare H 2 masers towards a NLS1 showing 
a low accretion rate and hence a smaller black 
hole mass is of great interest. According to a 
statistical study, it is found that the mean FIR- 
to-radio flux ratio of low-luminosity H 2 masers 
is higher than that of their high-luminosity coun- 
terparts, suggesting that the maser excitation of 
low-luminosity water masers such as the one in 
M82 is related primarily to star-forming activity, 
rather than AGN activity. 
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Table 1 

Summary of VLA Observations 



Object 


Date 


Tracking 


-Center 


Velocity 


AV 




Beam 


(T rms 


(Cal. Source) 


(d.m.y) 


R.A. 

(h m s) 


Decl. 

(° ' ") 


Range 
(km s" 1 ) 


(km s" 1 ) 


(hr) 


(HPBW) 

(mas) 


(mjy 
beam -1 ) 


M82 
(0955+654) 


13.4.02 


a 1950 

09 51 42.200 
09 54 57.847 


<5l950 

69 54 59.30 
65 48 15.53 


50-130 
70-150 


1.3 


9 


93x80 
PA=45° 


1.1 


M51 
(1419+542) 


4.7.03 


<*2000 

13 29 52.708 

14 19 46.597 


^2000 

47 11 42.79 
54 23 14.78 


365-525 
480-640 


5.3 


8 


82x77 
PA=6° 


1.1 


NGC4051 
(1153+493) 


4.7.03 


0:2000 
12 03 09.606 

11 53 24.466 


^2000 

44 31 52.52 
49 31 08.83 


640-725 
695-780 


1.3 


9 


77x72 
PA=81° 


1.3 



* Equinox for M82 is B1950.0; equinox for M51 and NGC4051 is J2000.0. 
a Velocity resolution 
b Total observing time 
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Table 2 
H 2 maser in M82 



Component 


OiB1950 


<5b1950 


Peak Flux 


Velocity 


Feature Width a 


Lh 2 o 




(09 h 51 m 00 B ) 


(69°54'00") 


(mJy) 


(km s~\LSR) 


(Channel) 


(Lq) 


Sl b 


42.6476 ± 0.0006 


58.180 ± 0.003 


13 


126 


2 


0.008 


S2 b 


42.1808 ± 0.0002 


59.263 ± 0.001 


03 


108 


3 


0.062 




42.1799 ± 0.0007 


59.257 ± 0.004 


15 


98 


1 


0.006 




42.1825 ± 0.0014 


59.252 ± 0.007 


8 


84 


1 


0.003 


Nl 


42.3951 ± 0.0012 


58.195 ± 0.007 


8 


105 


4 


0.007 


N2 


40.9436 ± 0.0007 


58.756 ± 0.004 


16 


146 


1 


0.006 


N3 


41.2053 ± 0.0014 


59.602 ± 0.007 


5 


98 


1 


0.002 



a Each channel is 1.3 km s 1 wide. 

b Adopting labels in lBaudrv fc Brouillet] U996D 
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Table 3 

H 2 Maser and Other Sources in M82 



H2O Maser 
(Position ID) 


Nearest Discrete Continuum a 


Origin a 
(09 h 51 m 00 B ) 


Nearest H92 c 


, b RRL 


CKS1950.0 

(09 h 51 m 00 B ) 


(5B1950.0 

(69°54'00") 


9 

(arcsec) 


OIB1950.0 

(69°54'00") 


(^S1950.0 


Sl(42. 64+58.1) 


42.694 


58.24 


0.31 


Hn 


42.61 


58.0 


S2(42. 18+59.2) 


42.210 


59.04 


0.22 


Hn 


42.19 


58.8 


Nl(42.39+58.2) 


42.481 


58.36 


0.49 


H11 






N2(40.94+58.7) 


40.938 


58.87 


0.17 


Hn 


40.97 


58.5 


N3(41. 20+59.6) 


41.302 


59.64 


0.51 


SNR 







IMcDonald et al.l l!2002fl 
' iRodriguez et al.l j2004D 

c Angular distance between each maser spot and the nearest continuum source. 
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Table 4 
H 2 Maser in M51 



Component 


a .72000.0 
(13 h 29 m 00 s ) 


£.72000. 

(47°11'00") 


Peak Flux 
(mjy) 


Velocity 
(km s~\LSR) 


Feature Width" 
(Channel) 


Lh 2 o 
(L ) 


Blue-shifted component 


52.5486 ± 0.0016 


43.337 ± 0.012 


8.2 


445 


2 


0.1 


Red-shifted components 


52.7089 ± 0.0001 


42.790 ± 0.001 


41.8 


539 


1 


0.47 




52.7087 ± 0.0001 


42.789 ± 0.001 


52.6 


565 


1 


0.58 




52.7091 ± 0.0001 


42.788 ± 0.002 


23.4 


576 


1 


0.26 



a Each channel is 5.3 km s 1 wide. 
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69 55 03 
02 
01 
00 
54 59 
58 
57 
56 

09 51 44 

69 55 03 
02 
01 
00 
54 59 
58 
57 
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09 51 44.5 



-i 
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I 






I 
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A 


S2 , 

+ * 

□ 

S1 + 




N3 

+ 


A 

N2 
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+ □ + a 
N1 


□ 

+ 




□ 








I 


I 






I 





44.0 



43.5 



43.0 42.5 42.0 41.5 

RIGHT ASCENSION (B1950) 



41.0 



40.5 



44.01+596 



44.0 



43.5 



43.0 



41.95+575 



42.5 42.0 41.5 
RIGHT ASCENSION (B1950) 



41.0 



40.5 



Fig. 1. — Top: Locations o f H2O masers, w ith other emission shown for comparison, in M82. Infra-red 
continuum peaks at 2.2 ilvci (Dietz et al. 19861 ) are plotted with stars, H2O masers are by crosses, 1667/ 1665 
MHz OH masers ( Argo et al. 2004 ) are shown by squares, 1720 MHz OH masers (Seaauist et al. 19971) are 



show n by triangles, and hard X-ray points are shown by broken crosses (jKaaret et al.ll2001uMatsumoto et al 
200 ll ). Larger squares and triangles indicate blue-shifted emission relative to 200 km s _1 ; their smaller 



counterparts show red-shifted (Vlsr> 200 km s _1 ) emission. The systemic velocity of the galaxy is 220 ± 5 
km s . It is obvious that all of the locations of OH and H2O masers are well separated in space, depending 
on their offset relative to the velocity of 200 km s -1 . Bottom: The 1.3-cm radio continuum version of the 
map shown in the top panel. Gray-scales are plotted from -1.0 to 5.0 mJy beam -1 . 
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Fig. 2. — Single-dish spectrum of the H2O maser towards the center of M82. The velocity resolution is 1.1 
km s _1 . The total integration time is ~ 60 minutes, observed with the MPIfR 100 m telescope at Effelsberg 
on 7 March 2002. The amplitude is scaled in janskys. 







Fig. 3. — VLA H2O maser spectra towards five locations in M82 within V t,br= 70 -150 km s -1 . The masers 
towards SI and S2 were originally detected at 1" resolution of the VLA ( Baudrv &: BrouilletJ 1996T ). Their 
positions are constrained at 0.1" resolution in this paper. The masers at Nl, N2, and N3 are new detections. 
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Fig. 4. — H 2 maser spectra of M51, observed with VLA-A on 4 July 2003. The left panel covers blue-shifted 
velocity range, and the right panel covers the red-shifted range. The velocity resolution is 5.3 km s _1 . The 
spectra were obtained from two IF channels, each with 12.5 MHz bandwidth. 
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Left: Large -scale 8.4 GHz nuclear continuum of M51 imaged by the VLA 



([Bradley. Kaiser, fc Baan II20041 ). superposed with two approximate positions of the H 2 maser. The red- 
shifted cluster is marked in red in the vicinity of the continuum peak and the blue-shifted feature(s) marked 
in blue is offset from both the peak and the other maser cluster. Right: First velocity moment map of 
the red-shifted maser from the uniform-weight map. A small velocity gradient over the maser is identified. 
Contours are spaced at 1-2 km s _1 and labeled with the absolute LSR velocity-scale. 
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Fig. 6. — H 2 maser spectra from the nucleus of NGC 4051 with a velocity resolution of 1.3 km s -1 using the 
VLA. The spectra consist of two independent IF channels that cover the velocity range of Vlsr= 640—780 
km s _1 , straddling the systemic velocity of Vlsr= 730km s _1 . 
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Fig. 7. — The 1.4 GHz radio flux density vs. FIR flux density (estimated from 60 /im and 100 /im). Crosses 
represent the known high-luminosity water masers in the literature, and filled circles represent low-luminosity 
wate r masers in star-forming nuclei (IC10 IC342 M51, M82, NGC253, NGC2146, and NGC69 46), primarily 
from lGreenhill et all (|l990l) . iTarchi et all (j2002al lbh ; iHenkel et all (|2002l) ; IClaussen et all (|l984l ) .Open circles 
represent prototypical OHMMs from Sanders et al. (1988). Mean constant ratios of FIR to 1.4 GHz radio 
flux densities for the three different samples are indicated by the solid line (high- luminosity H 2 maser), the 
dashed line (OHMM), and the dotted line (low-luminosity H 2 maser). 
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